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ABSTRACT 

The calibration of ground-based submillimetre observations has always been a difficult pro- 
cess. We discuss how to overcome the limitations imposed by the submillimetre atmosphere. 
Novel ways to improve line-of-sight opacity estimates are presented, resulting in tight rela- 
tions between opacities at different wavelengths. The submillimetre camera SCUBA, mounted 
on the JCMT, is the first large-scale submillimetre array, and as such is ideal for combating 
the effects of the atmosphere. For example, we find that the off-source pixels are crucial for 
removing sky-noise. Benefitting from several years of SCUBA operation, a database of deep 
SCUBA observations has been constructed to better understand the nature of sky-noise and the 
effects of the atmosphere on instrument sensitivity. This has revealed several results. Firstly, 
there is evidence for positive correlations between sky-noise and seeing and sky-noise and sky 
opacity. Furthermore, 850- (a.m and 450- (a.m sky-noise are clearly correlated, suggesting that 
450-|^m data may be used to correct 850-|xm observations for sky-noise. Perhaps most im- 
portant of all: if off-source bolometers are used for sky-noise removal, there is no correlation 
between instrument sensitivity and chop throw, for chop throws out to 180 arcsec. Understand- 
ing the effects of submillimetre seeing is also important, and we find that the JCMT beam is 
not significantly broadened by seeing, nor is there an obvious correlation between seeing and 
pointing excursions. 

Key words: submillimetre - instrumentation: detectors (SCUBA) - telescopes (James Clerk 
Maxwell Telescope) 



1 INTRODUCTION 

Ground-based observations at submillimetre wavelengths are 
severely hindered by the atmosphere, which absorbs, emits, and 
refracts the incoming radiation. High, dry sites are needed, such as 
Mauna Kea in Hawaii. Even then the transmission is generally poor, 
with only a small number of semi-transparent windows accessible. 

To calibrate submillimetre data, the atmospheric opacity must 
be accurately determined. The main absorber of radiation in this 
waveband is water vapour, although oxygen and ozone can be sig- 
nificant contributors. Furthermore, the transparency of the atmo- 
sphere often changes on short timescales. Thus, frequent measure- 
ments of the opacity are crucial, especially for the shorter wave- 
length windows at 350 and 450 \im. There is a strong dependence 
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of the transmission on wavelength; the shorter windows are more 
opaque, and deteriorate faster as conditions worsen. 

In addition to attenuating the signal, the atmosphere and im- 
mediate surroundings of the telescope emit thermal radiation sev- 
eral orders of magnitude larger than the source signal. Spatial and 
temporal variations in this sky emissivity give rise to "sky noise", 
which can degrade the effective instrument sensitivity by up to an 
order of magnitude. The thermal DC offset and sky noise variabil- 
ity can be largely removed by the conventional techniques of sky- 
chopping and telescope nodding. To be effective against sky noise 
requires the secondary mirror to switch between sky+source and 
sky faster than the rate at which the sky is varying (typically greater 
than a few Hz). Nodding the telescope primary to place the source 
alternately in both chop beams cancels slower varying sky gradi- 
ents due to chop-beam imbalances and time-dependent telescope 
spillover signals. It is not practical to nod the primary at the chop 



2 E. N. Archibald et al. 




400 600 800 1000 

Wavelength / ixm 

Figure 1. The SCUBA filter profiles measured in situ, superimposed on tlie submillimetre atmospheric transmission curve for Mauna Kea assuming 1 mm of 
precipitable water vapour. 



rate and so typical frequencies of 0.1 Hz are adopted. However, al- 
though these techniques diminish the effects of sky noise, they do 
not remove the residual signature completely. 

Fluctuations in the atmospheric refractive index, as the atmo- 
sphere drifts through the telescope beam, cause variation s in the 
path lengt h from source to telescope (|Church & Hills 199(]|; C hurch 
et al. 1993). This is sometimes referred to as submillimetre seeing, 
and causes apparent (short-term) pointing shifts, which ultimately 
degrade the S/N of an observation. This is particularly problem- 
atic for under-sampled arrays such as SCUBA. Until the time that 
adaptive optics becomes available for the submillimetre, the effects 
of this seeing must be closely monitored, and the data severely af- 
fected removed from a particular observation. 

This paper considers the limitations imposed on the quality 
of submillimetre data by the atmosphere, and presents methods 
that maximise the accuracy of the data. The layout is as follows: 
Section ^ provides a technical description of SCUBA and its ad- 
vantages compared with previous instruments, while Section ^dis- 
cusses techniques of determining atmospheric extinction as a func- 
tion of wavelength. We present our own skydip model and con- 
sider how continuously-operating radiometers are crucial for pre- 
cise measurements of atmospheric opacity. Sky-noise is described 
in detail in Section ^ particularly how it can degrade instrument 
sensitivity, and how most of its effects can be removed with state- 
of-the-art instrumentation such as SCUBA. The effects of refrac- 
tion noise (submillimetre seeing) are described in Section ^. Using 
a carefully constructed database of SCUBA observations and at- 
mospheric information. Section ^ analyses how the atmosphere di- 
rectly effects instrument sensitivity. In Section ^ we conclude with 
a discussion of future technology that will provide further improve- 
ments in the accuracy of submillimetre data. 



2 THE BENEFITS OF ARRAY RECEIVERS: SCUBA 

Until recently, observing in the submillimetre w as limited to single- 
pixel devices. UKT14 (Duncan et al. 1990), the forerunner to 
SCUBA on the JCMT, was detector-noise limited at all wavelengths 
of operation, except under periods of high sky variability, when the 
sensitivity at 350 and 450 \im was severely limited by sky-noise. 

An array receiver is a major improvement over a single-pixel 
system. The off-source pixels can be used to measure sky-noise on 
short timescales. This is crucial for accurate work: Omont et al. 
(1996) compared the performance of a detector array with a single- 
pixel system, and found that the poor sky cancellation offered by 
a single detector could in fact produce fake detections of faint 
sources. 

Within the past few years submillimetre astronomy has wit- 
nessed the arrival of the first large format bolometer arrays, together 
with detector sensitivities which are limited, under stable condi- 
tions, by the background photon noise from the sky. The largest 
and most powerful of this new generation of submm cameras is 
the Submillimetre Common-User Bolometer Array (H olland et al. 
1999), which operates on the 15-m James Clerk Maxwell Tele- 
scope, on Mauna Kea. SCUBA is a dual camera system containing 
91 pixels in the short- wavelength (SW) array and 37 pixels in the 
long-wavelength (LW) array. Background-limited performance is 
achieved by cooling the detectors to ~ 100 mK. Both arrays have 
approximately the same field-of-view (FOV) on the sky (2.3 ar- 
cminutes in diameter) and can be used simultaneously by means of 
a dichroic beamsplitter. The SW array is optimised for operation at 
450 \im, while the LW array is optimised for 850 \im. 

The wavelength of operation is selected by a bandpass filter 
carefully designed to match the transmission wi ndow. The fi lters 
are multi-layer, metal-mesh interference filters (Hazell 1991) lo- 



cated in a nine-position rotating drum that surrounds the arrays. 
They have excellent transmission (typically over 80%), and also 
less than 0.1% out-of-band power leakage. This latter characteristic 
is particularly important as it ensures that there is minimum con- 
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tribution to the source signal from extraneous sky emission. The 
spectral performance of the filters was measured by the Univer- 
sity rfLethbridge^^ourie^Jrans^^ et al. 
1994), and the resultant profiles are overlaid on the Mauna Kea at- 
mospheric transmission curve in Figure |lj For a more detailed look 
at the submillimetre atmosphere please refer to Vaylor et al. (2000 ). 



SCUBA skydips at 350, 450, 750 and 850 microns 

Skydip data from 7 Feb 1998 



SCUBA has three basic observing modes (H olland et al. 
1999). Photometry, uses a single bolometer to observe point 
sources. For sources larger than the beam, but smaller than the ar- 
ray FOV, a 64-point jiggle map produces a fully-sampled map for 
both the LW and the SW arrays, by jiggling the secondary mirror 
to fill in the undersampled arrays. Finally, sources larger than the 
array FOV are observed in scan-map mode. 

The SCUBA Upgrade Project was designed to improve the 
sensitivity of the instrument. The first part of the project was com- 
pleted in October 1999, and included the installation of two wide- 
band filters centred at 450 |i.m and 850 ^m (450W:850W). The 
wideband filters were designed to be more closely matched to the 
atmospheric windows (Figure and to be more sensitive than their 
narrowband predecessors (450N:850N) under all weather condi- 
tions. The measured improvement is a few percent at 850 ii.m and 
a factor of 2 at 450 (xm. The improvement at 850 (xm is not due to 
the width of the filter, but to the blocking filter that was installed to 
reduce contamination by infrared light. The overall spectral trans- 
mission characteristics of the 850-|i.m waveband are driven by the 
edge filters and detector feed-horn cut-off. Thus, unfortunately, the 
measured response of the 850- nm wideband filter is almost identi- 
cal to that measured for the narrowband filter. For the purposes of 
this paper, we will concentrate on the 450W:850W filters. 



3 ATMOSPHERIC ATTENUATION 

Determining the atmospheric attenuation of a source signal is crit- 
ical for calibrating submillimetre data. Assuming a plane-parallel 
atmosphere: 



(1) 



where Im and Jo are the signals incident at the telescope and at 
the top of the atmosphere respectively, A is the airmass (the se- 
cant of the zenith dista nce), and Tau is the zenith sky opacity 
( [Stevens & Robson 1994 ). Thus, precise measurements of Tau must 



be taken frequently. This is less crucial at 850 \xm; in good weather, 
Tau85o<0.3, and at low airmass, (<1.5), a 20% error in Taugso al- 
ters the measured source flux density by 5-10% at most. However, 
in worse conditions and particularly at 450 \im, an error in Tau can 
severely affect the measured source flux density. For example, as- 
suming a low airmass <1.5, a 20% error in Tau45o can alter the 
measured flux density by ~ 50 — 80%. 

Traditionally, Tau was derived by constructing a secant plot: 
the signal from the source is measured as a function of airmass, and 
assuming the sky does not change between measurements, the gra- 
dient of the plot gives Tau directly. However, this method requires 
a meaningful number of measurements of a bright source over a 
range of airmasses. For any level of acc uracy, the sky must remain 
very stable over a long period of time. Stevens & Robson (1994) 
demonstrate the difficulties of using this method to derive Tau. We 
describe here alternative methods that are able to track rapid varia- 
tions in Tau. 
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Tau(350)= 1.08+/- 0.012 




Figure 2. Typical skydip examples at 350, 450, 750 and 850 /im. Open cii- 
cles denote measurements taken with SCUBA; solid lines denote the model 
fits to the data. The skydips shown here correspond to a precipitable water 
vapour of ~ 0.7mm (see Equation 3). 



3.1 Skydip method 

SCUBA estimates the zenith sky opacity at the wavelength and az- 
imuth of observation by performing skydips. Skydips measure the 
sky brightness temperature as a function of elevation (usually be- 
tween 80 and 15 degrees), with absolute temperature calibration 
provided by hot and cold loads. The hot load is ambient temper- 
ature Eccosorb, and the cold load is a reflection of the cold op- 
tics inside the cryostat, with an effective temperature of ^ 60 K. 
An aperture plane chopper unit, spinning at 2 Hz, is used to switch 
rapidly between the sky and the two loads. The temperatures of the 
hot and cold loads are measured and corrected for the emissivity 
and reflectivity of the components in the optical path. A model de- 
scribing both the atmosphere (assuming a plane-parallel form) and 
the optical system is then fit to the data to calculate the zenith sky 
opacity: 



J rr 



= (1 — rjtel) Jtel + Vtel Jat 



bwf rjtel JatmS "^'"''^ (2) 



where Jmeas is the measured brightness temperature of the sky, 
ritei is the transmission of the telescope, Jtei is the brightness tem- 
perature of a black-body radiating at the temperature of the tele- 
scope, Jatm is the brightness temperature of the atmosphere, bwf 
is the bandwidth factor of the filter being used (1 — bwf is the 
fraction of the filter bandwidth that is opaque due to atmospheric 
absorption and, like Tau, is a function of water vapour content), Tau 
is the zenith sky optical depth and A is the airmass of the measure- 
ment. Technical details of this model are presented in Appendix ^ 
refer also to lHazen (199l[ ). 

Figure El shows typical skydips taken with SCUBA at several 
wavelengths. In practice, the skydip method provides an accurate 
measurement of Tau. However, it takes ~ 6 minutes to perform a 
skydip. Given this overhead, it is only practical to perform a skydip 
every 1.5-2 hours, and quite often the frequency is even less. An on- 
the-fly skydip mode has recently been commissioned, where data 
are taken contin uously, reducing the time taken to ~ 2 — 3 minutes 
dCoulson 200 lb. 
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where the pwv is in millimetres ( Davis et al. 1997 ) 
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Figure 3. CSO Tau data as a function of time (expressed as a fraction of 
the UT date) for two typical niglits. The CSO data are depicted by the '+' 
symbols. The 850-micron skydips taken for each night, s caled to 225 GHz 
values using the newly derived CSO Tau relations(Section [3.3| ), are denoted 
by the solid circles. The solid line is a polynomial fit to the CSO Tau data. 
The fit residuals are also shown. 



3.2 CSO Tau Monitor 

The nearby Caltech Submillimetre Observatory (CSO) operates a 
225 GHz (1.25 mm) tipping radiometer, which performs a skydip 
every 10 minutes, albeit at a fixed azimuth. The precipitable water 
vapour of the atmosphere is related to the CSO Tau as follows : 



Comparing skydips taken with SCUBA to those taken at the 
CSO yields relations between Taucso & Taugso and Taucso & 
Tau45o. If the scatter about these relations is small, the CSO Tau 
monitor can be used to measure the opacity more frequently than 
we can with SCUBA, with no additional overhead. 

The CSO data show a significant amount of noise, but track 
long time-scale (~2 hour) variations very well. This noise is more 
than one would expect given the measurement errors quoted in the 
CSO Tau archive. It is probable that this is due to instrumental noise 
and does not represent the behaviour of the sky. A polynomial can 
be fit to the data to track the large-scale variations in CSO Tau as 
opposed to the (presumably) instrument noise. 

In Figure|3|, example polynomial fits are presented for two typ- 
ical nights. Using the relations derived in Section 3.3, we have also 



pwv = 20{Taucso - 0.016) 



(3) 



plotted the 225 GHz Tau predicted by the 850-tim skydips taken 
each night. It is striking how well the skydips track the polynomial 
fits, even when one would imagine that the CSO Tau was moving 
around too much to be useful. 

Producing a composite picture of the night, using the CSO 
Tau, the SCUBA skydips, and the polynomial fit, provides an addi- 
tional level of quality control. It can give the observer a much bet- 
ter feel for how the atmosphere was actually behaving on a given 
night, especially if conditions were apparently unstable. Consider 
first the lower plot in Figure ^. It is clear that interpolating between 
the SCUBA skydips would give an erroneous measurement of Tau. 
In addition, the extreme scatter in CSO Tau indicates that the end 
of the night was unusable. 

The composite picture can provide further information about 
when data should be treated with care. For example, on some nights 
the CSO Tau and SCUBA skydip data disagree with each other, on 
other nights both display unusually high levels of scatter, indicat- 
ing an inherently unstable night. If, on the other hand, the CSO 
Tau has a high-level of scatter but the SCUBA skydips follow the 
polynomial fit, the scatter is unlikely to be representative of the sky 
itself. 



3.3 Tau relations 

The time-resolution of the polynomial fits to the CSO Tau is im- 
pressive, ~ 2 minutes if the residuals are small, as they typically 
are. We can capitalise on this by deriving tight relations between 
the CSO Tau and the opacity measured at SCUBA wavelengths. 

These relations have been constructed taking the following 
points into account: 

(i) Skydips were discarded where the model failed owing to, for 
example, the atmosphere losing its plane-parallel nature or a cloud 
drifting overhead. In these cases the model returns unrealistic val- 
ues for the fit parameters and/or unusually high fit residuals. This is 
a more common occurrence at 450 [im, where the shape and height 
of the atmospheric window are more sensitive to the presence of 
water vapour (see Figure |^). We discarded ~20% of the data at 
850 ^tm and ~50% of the data at 450 vim. Note, we can still cali- 
brate at 450 vim when the skydip fit fails as we can extrapolate from 
either the CSO or the 850- vim Tau. 

(ii) We ignored data taken when the CSO Tau monitor was off- 
line. 

(iii) If either the CSO Tau monitor or the skydip indicated a non- 
physical value of Tau, i.e. negative or zero, the observation was 
ignored. 
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Figure 4. The wideband Taugso-Taucso relation. For the left plot, the CSO Tau data were taken directly from the archive, and every skydip observation was 
accepted. The right plot shows the relation with Taucso calculated from the polynomial fits, and with the poor-fitting skydips having been removed. 
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Filter System 


Time Period 


Tauy 


Taux 


a 


b 


450N:850N 


Feb. 04, 1998-Oct. 10, 1999 
(pre-upgrade) 


Taugso 
Tau45o 
Tau45o 


Taucso 
Taucso 
Taugso 


3.99±0.02 
23.5±0.2 
5.92±0.04 


0.004±0.001 
0.012±0.001 
0.032±0.002 


450W:850W 


Dec. 05, 1999-Sept. 30, 2000 
(post-upgrade) 


Tau850 
Tau45o 
Tau45o 


Taucso 
Taucso 
Tau85o 


4.02±0.03 
26.2±0.3 
6.52±0.08 


0.001 ±0.001 
0.014±0.001 
0.049±0.004 



Table 1. Tau relations for both the narrow and wide 450:850 filter systems, calculated using least-squares regression. The corresponding errors are 1 — ct. The 
relations have been constructed using the CSO polynomial fits and by ignoring poor-fitting skydips. 

(iv) The datasets were restricted to Taucso < 0.2 (SCUBA is 
not used in weather conditions worse than this). 

(v) A model of the form Tauy ~ a{Taux — b) was fit to the 
data to derive the Tau relations. 

Figure ^ shows how much the scatter in the relations has been 
reduced by simply using the polynomial fits to estimate the CSO 
Tau and ignoring untrustworthy skydips. The final relations are pre- 
sented in Table |l| and Figure ^. They display relatively little scatter, 
even in the poorest weather conditions, and can thus provide accu- 
rate submillimetre calibration. 

Comparing the wideband and narrowband filters, the wide- 
band CSO Tau relations are steeper at 450 vim but are almost iden- 
tical at 850 p.m. At first sight, the difference at 450 p.m is perhaps 
unexpected given the lower central wavelength of the wideband fil- 
ter (at 850 vim, the narrow and wideband filters have almost identi- 
cal central wavelengths). However, the 450 |i.m wideband filter in- 
cludes an H2O line that is always present, even in very dry weather 
conditions. Thus, the wideband filter will always have a steeper 
slope than the narrowband filter. 

If the submillimetre opacity is due solely to water vapour, 
the Tau relations for the different filters are expected to intercept 
the origin: if there is no water vapour in the atmosphere, Taucso, 
Taugso and Tau45o will all equal zero. However, we have found ev- 
idence for non-zero intercepts. Assuming the straight-line model 
can be extrapolated to the intercept, this could be explained as fol- 
lows: the relative contribution of ozone vs. water vapour to absorp- 
tion is greater at 225 GHz than at 850 vim and in turn greater at 



850 vim than at 450 vim. The ozone contribution is relatively invari- 
ant for long periods of time, and these offsets should be constant. 



4 SKY-NOISE 

4.1 Signatures of sky-noise 

The atmosphere emits submillimetre radiation several orders of 
magnitude larger than the signals we are trying to measure. Varia- 
tions in this emission result in sky-noise, which manifests itself as a 
1/f component in the noise spectrum. This is illustrated in Figure]^ 
where the 1/f component extends out to 4 Hz and is well above the 
system noise level. Two-position chopping at higher frequencies, 
~ 8 Hz, greatly reduces this sky-noise contamination (D uncan et al. 
I995I 

As noted earlier, telescope nodding is also essential for the 
removal of sky-noise. This can be illustrated by the Allan variance, 
which is simply the variance of the difference of two contiguous 



measurements (e.g., Allan 1966; Schieder & Kramer 2001). When 



plotted, a slope of —1 indicates white noise, a slope of denotes 
1/f noise, and a slope of -l-l is drift noise. Figure ^ shows the Allan 
variance data as a function of time for a typical observation with 
and without nodding (in both cases the telescope was chopping at 
~ 8 Hz). A clear 1/f component is present if the telescope is not 
nodded: after around 20 seconds the noise stops integrating down. 
This component disappears when nodding is employed. 

However, it is impractical to nod the telescope at the chop 
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450Tau=6.52(850Tau-0.049) 
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Figure 5. Tau relations for the wideband 450W:850W filter system. The top 
plot depicts the relationship between the SCUBA skydips and CSO Tau. The 
bottom plot depicts the Tau45o-Taug5o correlation derived by comparing 
450- |jm and 850- [am skydips. Models of the form Tauy = a(Taux — b) 
have been fit to the data in every case. 



frequency. Thus, even with a chop/nod configuration, short-term 
temporal variations in the sky emissivity, as well as spatial effects 
caused by chopping through slightly different atmospheric paths, 
still exist. This residual sky-noise must be removed. There are sev- 
eral ways to do this; filters can be designed to select the most trans- 
parent parts of the atmospheric transmission window. The optical 
throughput should be single-moded, i.e. the minimum required to 
couple to a point so urce with maximum spatial resolution and min- 
imum background (Duncan et al. 1995). Keeping the chop throw 
as small as possible, and in an azimuthal direction (i.e. chopping 
through the same atmospheric layer), should also help. In this pa- 
per we will concentrate on the benefits of having a bolometer array 
to aid the removal of sky-noise. 
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Figure 6. Noise spectrum for LWA central pixel, showing 1/f noise tail of 
sky-noise. The system noise trace is from the central bolometer looking at 
a 4 K blank placed over the arrays. 
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Figure 7. Allan variance plot showing 1/f noise, white-noise and drift noise. 
4.2 Removal of sky-noise 

In the submillimetre, the emissivity variations o ccur in atmospheric 



cells that are large r than the array field-of-view ( lenness et al. 1998 



Borys et al. 19%). Therefore, the noise should be correlated to a 



large extent across the array and also between wavelengths. 

One of the main advantages of SCUBA is that it contains an 
array of bolometers. When observing a compact source, the off- 
source bolometers can be used to measure the sky signal on short 
timescales and hence remove any residual sky-noise that chopping 
and nodding have failed to account for. 

This is clearly illustrated in Figure ^ where the signal of both 
the centre on-source 850-n.m pixel and the sky (measured by av- 
eraging the outer ring of bolometers) are shown. The signals are 
highly correlated, and subtracting the sky from the source reveals 
a clear positive signal (in which the jiggle modulation is appar- 
ent). Using this sky-subtraction method can improve the S/N sig- 
nificantly. In Figure ^ we also display the 450- jtm centre pixel, and, 
as expected, it is also highly correlated with the 850- [xm data. 
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Figure 8. Example of the identification and removal of sky-noise. This fig- 
ure first appeared in Holland et al. (1999) - note one integration comprises 
128 seconds, split into 4 exposures of 16 jiggle points, with 1 second ob- 
served per jiggle point per nod position (there are two nod positions). 

Failure to remove this residual sky-noise can, in the worst 
conditions, degrade the sensitivity of the instrument by an order 
of magnitude. Figure ^ displays the noise integrating down with 
time for a deep 850- (am photometry observation. Data are shown 
for the standard chop/nod configuration, and also after the residual 
sky-noise has been removed using the off-source bolometers. We 
define the sensitivity, also known as the noise equivalent flux den- 
sity (NEFD), as the noise reached in 1 sec of integration. It is clear 
that the NEFD of the instrument is significantly better if the residual 
sky-noise is removed. Furthermore, at the end of the observation, 
there was a very obvious, large deterioration in the sky conditions. 
This is plainly visible in the case where the residual sky-noise was 
not removed. Removing the residual takes this change of conditions 
into account, and allows the noise to integrate down oc ^/t. 

For scan-map observations, removing the sky-noise is more 
difficult, as every bolometer could be observing either the source or 
the sky at a given time. The source can be removed from the datas- 
tream by the simple assumption that the source structure is constant 
over time, while the sky is varying on timescales of a few seconds. 
The sky-emission noise can then be calculated and removed from 
the data (Jessop et al. in preparation). 

4.3 Investigation of sky-noise 

For the purposes of this paper, we define sky-noise as the standard 
deviation of the sky signal, with any errors associated with nodding 
removed. The dataset used here consists of deep photometry obser- 
vations only, as it is easy to specify the sky bolometers and the long 
observations provide good statistics. 

The relationship between sky-noise, seeing, and sky opacity 
is not immediately obvious; the scatter in the data was too large to 
reveal an obvious trend. To overcome this, we binned the data. This 
also posed problems as the data were not Gaussian distributed and 
it was not clear whether the mean (which is skewed by outliers) 
and associated standard deviation (which is unnaturally large given 
the scatter in the data) were useful quantities to measure. Instead 
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Figure 9. Standard error evolution with time for a deep 850- (am photometry 
observation. For the open circles, only chopping and nodding have been 
used for sky cancellation. For the solid circles, the sky-noise residual has 
been removed using the off-source bolometers. This data was, for the most 
part, taken in uncommonly stable weather conditions. 

we measured the mode and full-width half-maximum of the data 
in each bin. These plots are shown in Figure |lo| and although the 
'errors' are somewhat large, there is evidence for positive corre- 
lations between sky-noise and seeing and between sky-noise and 
Tau. There is, however, no evidence to support the long-standing 
belief (from single-pixel instruments) that sky-noise often increases 
in very dry weather conditions. 

Considering the 850- nm and 450- urn sky signals for 20 
minute photometry observations. Figure hll presents a histogram 
of the Spearman-rank correlation coefficients. A value of -l-l indi- 
cates a pure positive correlation, indicates no correlation, and - 1 
indicates a pure negative correlation. On the whole, the 850- nm 
and 450-^im data appear to be highly correlated. However, there 
are clearly times when the strength of the correlation is low. 

Given the known correlation between sky-noise at 850 [im and 
450 ^Lm, it has been suggested that the short wavelength array could 
be used to provide sky cancellation at 850 \Lm. This would be ex- 
tremely useful for faint sources where the 450 [im field-of-view is 
likely to be source-free, but the 850 (xm field is expected to contain 
several sources or extended structure. 

The ratio of the 450- [xm sky signal to the 850-|J.m sky signal 
(corrected for extinction) is displayed in Figure |^. For the narrow- 
band data, the median ratio is ~ 4.5 in agreement with Borys et al. 
(1999) and [Tenness et al. (1998| ); for wideband data, the median 
ratio is ~ 8 and the distribution is wider. The ratio shows less dis- 
persion if it is extinction corrected than if not, suggesting that the 
water vapour creating the sky-noise is higher up in the atmosphere 
rather than just above the telescope (where zero correction would 
be more suitable). 

In Figure |l^ we present the ratio of 450-|i.m:850-|i.m sky- 
noise as a function of the corresponding transmission ratio. Thus, 
for a given transmission ratio, it is possible to use this plot to cor- 
rect 850-|J.m data for sky-noise using 450- ^im sky-noise data. It is 
clear that the sky-noise is larger for the wideband filter set than for 
the narrowband. This difference is to be expected as the 450- vim 
wideband filter is double the bandwidth of its narrowband prede- 
cessor, and contains an H2O absorption line, while the two 850- 
li.m filters have almost identical bandwidths (Duncan 1983 showed 
analytically that for a sky-limited system NEP oc Ai^). This ex- 
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Figure 10. 850- mn .sky-noise binned in CSO Tau (top plot) and seeing (bot- 
tom plot). The data points are the mode of each bin, and the 'errors' are the 
corresponding full-width half-maxima. 



— All Data: 72.1 % have S-R > 0.75 

-- Wideband Data: 75.9 % have S-R > 0.75 

Narrowband Data: 68.3 % have S-R > 0.75 
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Figure 12. Histogram of the ratio of 450- (Jin sky signal to 850- [om sky 
signal. 
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Figure 13. The 450- |am:850-|am sky-noise (not corrected for extinction) ra- 
tio against the 450- |am:850-|am transmission ratio. The black circles denote 
the narrowband filter, the gray circles denote the wideband filter 



cess sky-noise is the cost of the improvement in sensitivity offered 
by increasing the 450-|i.m bandwidth. However, as we have already 
demonstrated, this excess is removable with a bolometer array. 

A transmission ratio of 1 in Figure |n| corresponds to there 
being no atmosphere between the telescope and the source of the 
sky-noise. If the sky-noise ratios in Figure [l^ are extrapolated to 
a transmission ratio of 1, we find values of ~ 3.5 and ~ 6.5 for 
the narrowband and wideband filters correspondingly. These num- 
bers are ingood agreement with the theoretical values derived in 
Appendix ^ The extinction-corrected ratios in Figure are larger 
however, indicating that we are over-correcting for extinction in 
calculating these ratios. This in turn suggests that the water vapour 
responsible for the sky-noise lies somewhere in the middle-upper 
layers of the atmosphere, not at the very top of the atmosphere it- 
self. 



an Rank Correlation Coeffioient 



Figure 11. Histogram of Spearman-rank correlation coefficients comparing 
850 yim and 450 |am sky signals. 



5 SUBMILLIMETRE SEEING 
5.1 Measurement on Mauna Kea 

Submillimetre seeing arises from variations in the refractive index 
of the atmosphere, principally due to the passage of water vapour 
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Figure 14. Fractional signal as a function of pointing drift, assuming a 
FWHM of 14.2" at 850 |am and 7" at 450 [am. 



through the beam. This variation is measured on Mauna Kea by a 
phase monitor operated by the Smithsonian Astrophysical Obser- 
vatory (SAO). The SAO device comprises two 1.8 m dishes placed 
about 100 m apart located near the JCMT The dishes point low in 
the eastern sky, and detect a signal at 12 GHz from a geostationary 
satellite. The difference in path lengths results in a phase differ- 
ence between the two received signals that changes slowly with the 
oscillation of the satellite position. Turbulence in the atmosphere 
adds noise to the phase difference. The phase difference is mea- 
sured each second and one minute's worth of data are analysed for 
the rms scatter to provide the seeing at 12 GHz over a 100 m base- 
line. 

The conversion from the rms of these phas e fluctuations to an 
rms seeing for the 15 m JCMT is suggested by Masson (1991 ) as 
being: 



JCMT seeing (arcsec rms) — 0.5 SAO seeing (degrees rms) 

Fifteen minute averages of these seeing data are available at the 
JCM T as part of the Telescope Management System (T ilanus et al. 
1997). By the middle of the night the JCMT seeing typically drops 
to 0.2 arcsec, and observers consider seeing > 1 arcsec to be 'poor'. 



5.2 Impact of seeing: broadening & displacement 

If atmospheric turbulence on scales smaller than the telescope di- 
ameter is signi ficant then beam-broade ning will occur. However, 
experiments by Church & Hills (1990) indicate that these effects 



are small at the JCMT. Beam displacement results from large- 
scale, anomalous refraction, or 'tip/tilt' effects, and effective beam- 
broadening results if observations take longer than the timescale of 
these fluctuations; the time average of the beam motion being su- 
perimposed on t he beam pro flle. Anomalous refraction timescales 
are of order 1 s (Olmi 2001), which is considerably shorter than 



standard JCMT observations (minimum 18 s), but comparable with 
the SAO monitor rate. JCMT jiggle maps (integration times > 32 s) 
therefore suffer both minimal average pointing shifts and minimal 
image broadening as the turblence tends to occur on scales larger 
than the telescope diameter. 

Additionally, analysis of JCMT tracking data shows that there 
is no correlation between seeing and pointing excursions. Tracking 
data are long jiggle maps of bright point sources, and each integra- 
tion (32 s worth) is analysed for the location of the image centroid 



to yield information on pointing offsets as a function of azimuth, 
for instance. Some tracking datasets cover a sufficiently long time 
interval that significant changes in seeing occur. Within such data 
there is no correlation between the change in seeing and excursions 
in pointing. 

Tracking data also tend to validate the conversion factor (0.5) 
between the SAO phase (rms, degrees) and JCMT-seeing (rms, arc- 
seconds), since the scatter of the pointing offsets about some low- 
order polynomial fit is usually of the same size as the JCMT-seeing. 



5.3 Do errors in computed refraction generate pointing 
errors? 

Relatively thin atmospheric layers (thicknesses < 100 m, say) with 
arbitrary values of temperature, pressure, and humidity have no im- 
pact upon the overall refraction: the refraction towards the normal 
of a ray entering such a layer is nullified by the refraction away 
from the normal upon exit. Image displacement (a change in point- 
ing) must therefore result from either non-linear atmospheric struc- 
tures or from uncompensated changes in local humidity. 

The latter hypothesis is simple enough to test using archived 
JCMT/SCUBA pointing data, and no such relationship is found. 
The formulae in Appendix ^ show that, to first order and at con- 
stant zenith distance, elevation pointing corrections of 0.0768(/i — 
20) arcseconds need to be applied to account for the impact of lo- 
cal humidity, h(%). In the absence of this correction a plot of the 
change in elevation pointing against the change in humidity would 
be expected to have a slope of 0.0768, whereas a plot of the same 
data perfectly corrected for humidity will have a slope of zero. 
Analysis of 3000 pairs of consecutive pointing measures yields a 
slope of —0.009 ± 0.018, commensurate with the latter scenario. 
The source of residual pointing shifts therefore would seem to be 
excursions from the idealized atmospheric model - i.e. turbulence - 
which is probably no great surprise. 



5.4 Impact of pointing errors upon flux measurements 

Whatever their origin, pointing errors of size 9 (arcseconds) cause 
a signal of strength So to be measured as 



S = So exp 



-4 ln{2) 



FWHM 



(4) 



where FWHM is the full width at half-maximum of the Gaussian 
beam. These are shown graphically in Figure |l4| in the cases of the 
JCMT operating at 850 |i.m and 450 ^.m. 

As to the distribution of pointing errors as a result of refraction 
noise, we have determined that the expected signal, < 5 >, is 
related to the true peak signal. So, by 

So/ < S > = 1 + 0.0275 for a 14.2" beam at 850 fj,m 

So/ < S > =1 + 0.1132 for a 7" beam at 450 

where a is the rms refraction noise. The full derivation of these 



equations can be found in Appendix C2 



6 THE EFFECT OF THE ATMOSPHERE ON 
SENSITIVITY 

We have recently made a concerted effort to characterise SCUBA's 
performance. One aspect of this work was the creation of a database 
of calibration observations. The key result from studying this data 
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Figure 15. SCUBA sensitivity as a function of normalised sky transmis- 
sion for tlie 850W (black circles), 850N (open blue circles), 450W (green 
squai'es), and 450N (open pink squares) filters. For the top plot, only chop- 
ping and nodding were used for sky cancellation. For the bottom plot, the 
off-source bolometers were also used for sky cancellation. 



is that the flux conversion factors (FCFs) ar e quite stable, to a n 
accuracy of 5% at 850 \im and 20% at 450 vim ( Jenness et al. 2001 ). 
Using this database, we have created a homogeneous dataset of ~ 
4000 calibrated photometry observations, to better understand the 
effect of the atmosphere on sensitivity. 

Figure |l5| presents the NEFD as a function of normalised sky 
transmission for a chop/nod configuration, and also using the off- 
source bolometers to remove the residual sky-noise. The NEFD is 
significantly lower and is inherently more stable if the off-source 
bolometers are used for sky-cancellation. Data are shown both for 
the wideband filters, and the pre-upgrade narrowband filters. It 
is worth noting that the 450-nm wideband filter (which is better 
matched to the entire atmospheric window but includes an H2O 
line) is considerably more sensitive under all conditions. There is 
little difference between the narrow and wideband 850-|i.m filters, 
which is to be expected because, as noted in Section ^, the over- 
all filter transmission profiles are almost identical as measured in 
SCUBA. 

We have also considered how the sensitivity of the instru- 
ment is affected by the chop throw. Figure presents the data 
for two cases: a standard chop/nod method of sky cancellation, and 
a chop/nod method where the off-source bolometers are used to 
remove residual sky noise. Although the error bars are somewhat 
large, the standard chop/nod method seems to support the long- 
held belief that sensitivity decreases with chop throw (e.g. Church 
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Figure 16. NEFD vs. chop throw for both the wide and narrowband 850- 
[om filters. The dataset is restricted to the normalised sky transmission lying 
between 0.5 and 0.75, and the values are normalised to a transmission value 
of 0.65. The top plot is for a chop/nod method of sky cancellation and is 
restricted to observations with an azimuthal chop throw. In the bottom plot 
off-source bolometers ai'e used to remove the residual sky-noise, and there 
is no restriction on the chop throw angle. 



et al. 1993; Duncan et al. 1995). However, when using the array 
and when off-source bolometers are used for sky cancellation, large 
chop throws do not degrade the sensitivity. In this case, the NEFD 
vs. chop throw plot is perfectly flat out to a chop throw of 180 arc- 
sec. The only variation is due to an increase of the FCF with chop 
throw, of the order of 10% between 45 and 120 arcsec (J enness 
etal. 2001). 



7 SUMMARY AND THE FUTURE 

As described in this paper, we have a good grasp of the atmospheric 
limitations in the submillimetre and how best to overcome them. 
Our main findings can be summarized as follows: 

(i) Tau relations between the CSO Tau at 225 GHz and both the 
850 ii.m and the 450 \im. SCUBA filters. The relations display rela- 
tively little scatter. 

(ii) For accurate sky cancellation, it is essential to use off-source 
bolometers to monitor and remove the sky-noise. 

(iii) There is evidence for positive correlations between sky- 
noise and seeing, and sky-noise and sky opacity.but the correlation 
coefficient varies significantly depending on the dataset Figure |ll]. 

(iv) 850- (i.m and 450- fxm sky-noise are clearly correlated, but 
there are times when the correlation is low. 
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(v) The JCMT beam is not significantly broadened by seeing. 

(vi) There is no obvious correlation between seeing and pointing 
excursions. 

(vii) If off-source bolometers are used for sky cancellation, 
chopping as far as 180 arcsec (in any direction) does not affect the 
sensitivity of the instrument. 

In the future, as our understanding grows, as technological ad- 
vances are made, and as the next generation of submillimetre tele- 
scopes are destined to higher and drier sites, we will be able to do 
even better. Some things to look forward to are: 

(i) Line-of-sight radiometers which give real time estimates of 
Tau in the direction you are observing. One of these radiometers is 
now in operation at the JCMT (e.g. Wiedner et al. 2001, Phillips et 
al. in preparation). 

(ii) Instantaneous estimates of the line-of-sight opacity should 
be possible using a hot load, cold load, and sky observation at the 
elevation of the source, given a database of sky temperature versus 
elevation information (e.g. Smith, Naylor & Feldman 2001 ). 

(iii) High-speed data sampling leading to a more efficient ob- 
servingjnode_called_DREAM|_^Le Poole & van Someren Greve 
1998). In DREAM mode, the functions of chopping and jiggling 
the secondary mirror are combined to a single step action, eliminat- 
ing the need to sample empty sky for half the time. Another gain 
from high-speed sampling is real-time suppression of sky-noise. 

(iv) Adaptive optics for the submillimetre. The Large Mil- 
limetre Telescope (LMT) is currently designing a radiometric 
wave-front sensor which will measure the tilt of the incom- 
ing wavefront to compensate for atmosphere-induced pointing er- 
rors (For further de tails refer to the LMT website: http://www- 
Imt.phast.umass .edu/) . 

(v) DC-coupl ed fully sampled arrays such as SCUBA 2 (H ol- 
land et al. 2000; [Robson et al. 2001 ), which would remove the 
necessity to chop and nod for sky-noise removal. This is an ad- 
vantage: there is no chance of chopping onto a nearby source, and 
there will be a sensitivity gain as all the time is spent looking at 
the source instead of half the time looking at the sky. In addition, 
chopping limits the size-scales we see in maps to be no more than a 
few times the chop throw. Finally, image reconstruction techniques 
for two-beam chopping tend to propagate noise, and so with only a 
single-beam on the sky this problem will be minimsed. 

(vi) In the future, submillimetre telescopes will be built at sites 
where the atmosphere is largely transparent, for example Chajnan- 
tor in Chile at 17,000 ft. (ALMA) and the South Pole. The extreme 
cold at the South Pole results in very small amounts of water vapour 
i n the atmos phere, and thus sky-noise is much lower than other sites 
(^ark 200 1|). 



(vii) Interferometers such as ALMA have a unique advantage: 
the atmospheric signals decorrelate, and therefore sky-noise should 
not be a problem. 
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APPENDIX A: SKYDIP MODEL 

When performing a skydip, SCUBA measures the sky brightness 
temperature as a function of airmass. W e present a multi-layer 
model of the atmosphere (Hazell 1991) which, when compared 
with the data, yields the zenith sky opacity. The model at each 
wavelength takes the form: 



J rr 



— (1 " Vtel) Jtel + Vtel J atn 



' bwf ritel Jatme ^''"'^(Al) 



where Jmeas is the measured brightness temperature of the sky, 
rjtei is the transmission of the telescope, Jtei is the brightness tem- 
perature of a black-body radiating at the temperature of the tele- 
scope, Jatm is the brightness temperature of the atmosphere, hwf 
is the bandwidth factor of the filter being used (1 — bwf is the 
fraction of the filter bandwidth that is opaque due to atmospheric 
absorption and, like Tau, is a function of water vapour content), Tau 
is the zenith sky optical depth and A is the airmass of the measure- 
ment. 

Of these parameters, Jmeas, Jtei and A are known. Jatm can 
be estimated from the ambient air temperature at ground level us- 
ing a model for the behaviour of the observing layer above the tele- 
scope, as described below, rjtei may be fitted to the data for every 
skydip and, because it does not vary with atmospheric conditions, 
a reliable 'average' value can be derived from many observations. 
Thus, there are two remaining free parameters, r and hwf, that 
must be derived from the fit. 

Jatm is calculated from Tamb, the ambient air temperature, 
by assuming that the sky emission is dominated by a single ab- 
sorber/emitter whose density falls exponentially and temperature 
linearly with height. In this case it can be shown that: 

Jatm = Jamb Ak exp ^ Ji2^ ^ 

exp [Akh2 (exp (-i^) - l)] (^1 - ^) dh (A2) 

where hi is Jamb/()-5 to give a 6.5 K fall in temperature per km 
height, /i2 is the scale height of the absorbers (2 km), A is the air- 
mass and k the extinction per km. 

If we approximate the result of the integral by: 



Jatm = JambXg [1 - GXp {- Akh'i, 



^9 = 1+ rj. exp 

J- amb 



ATau 

qconst 



(A3) 



(A4) 



where Tiapse is the temperature drop per kilometre altitude 
(—6.5 K/km) and Xgconat is a constant determined empirically and 
has a value of 3.669383. 



APPENDIX B: THEORETICAL DERIVATION OF 
SKY-NOISE RATIO 

If we consider two lines of sight separated by an angle theta and 
identical to each other except for the fact that one line of sight has 
a cloud of water vapour at a distance to which the line of sight 
opacity is Tdoud, and the cloud itself has optical depth 5t, then 
the difference in surface brightness between the two lines of sight 
(measured in units of brightness temperature) is simply 



5 J = eXp{ — Tcloud) St {Jcloud — Jinc) 



(Bl) 



where Jdoud is the physical temperature of the cloud, and 
Jinc is the incident flux on the cloud from the upper atmosphere. 
In all cases it is likely that Jdoud is greater than Jinc so that the 
observed effect of the cloud is to increase the surface brightness 
along that line of sight. 

If the cloud is close to the top of the water vapour column den- 
sity then Jinc is negligible compared to Jdoud and one can simply 
write: 



5 J = exp{— Tdoud) St Jdoud 



(B2) 



If one considers the signal at two wavelengths Al, A2 then the 
ratio is 



SJxi _ exp{-Txi) Stxi 
5J\2 exp(-Tx2) Stx2 ' 



(B3) 



where we have removed the subscripts cloud. One can express 
this in terms of power by simply invoking the Rayleigh- Jeans law 

B{X) cx J/{X^) so that, 



SBxi _ exp{-Txi) Stxi 



SBx2 exp{-TX2) Stx2 A22 ' 

and also to measured volts by applying an appropriate FCF 

SVx^ ^ exp{-Txi) Stxi Al" FCFx2 
SVx: 



exp{-Tx2) Stx2 A22 FCFxi 



(B4) 



(B5) 



In the case where the opacities are simply proportional to wa- 
ter vapour column density (as we have shown in this paper) then 
Stxi/Stx2 is a constant, and we expect the ratios of the signal at 
the two wavelengths to be proportional to the ratios of transmission 
from the cloud to the telescope at the two wavelengths. One also 
expects 



SJx 



SJx 



(B6) 



exp[-Txi) exp{-Tx2) 

where one has to be careful to note that the value of r is to the 
cloud, not all the way to the top of the atmosphere. 

For the SCUBA filters the ratios of 5r850ri/(5T450n and 
Stssoiv /St 4sow have been measured earlier in this paper. J enness 
et al. (2001) have measured the FCFs from point sources using a 
60arcsecond radius aperture, which are also (by symmetry) valid 
for a uniform extended source of size 60 arcseconds. In the section 
on sky noise in this paper most of the photometry observations were 
taken with 60 arcsecond chop throws so that the skynoise measured 
is a measure of the difference in signal between two largely (but not 
exactly) overlapping regions of diameter 15 meters, 60 arcseconds 
apart. If the cloud is at a height, h above the telescope then the di- 
ameter of the beam in steradians is 15/h(m) or 25 arcminutes at a 
height of 2 km. This explains why the sky noise is uniform over the 
array. 

In Fourier space we can consider the beam as a high spatial 
frequency filter of the form 
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exp 



15 



(B7) 



i.e. one that filters out spatial scales < 15 m and the chop as a low 
spatial frequency filter of the form 



sin(2 TT k ch) 



(B8) 



where k is the wavenumber, c is the chop throw in steradians, and 
h is the height of the cloud in metres. This can be simplified to 
2nkch for the frequencies that are passed by the beam. The effect 
of having two differing spatial frequency filters is to couple only 
weakly (at a level ch/15) to clouds of approximate spatial scale 
15 m. Clouds smaller than this are filtered out by the beam, clouds 
larger than this are filtered out by the chop throw. For typical values 
(h = 2 km, c = 60 arcseconds), the coupling is only 4 percent - 
illustrating how effective chopping is as a first stage to removing 
sky noise. 

Althougl^_strictl^_s£eaking_^he_^CF^ 

et al. (2001) are only valid for extended far field objects, one can 
use the ratio of the values measured as an estimate of the ratio of 
the FCFs which would be valid for nearfield extended objects (i.e. 
the cloud in question). The validity of this assumption can then be 
assessed by comparison with the data. 

If one therefore uses the various values discussed in the previ- 
ous paragraph then one finds: 



SVisow 

SVs50W 

SV450N 
and 

SBssow 
SB450N 



6.64 



3.32 



: 23.2 



21.2 
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-T450w) 


exp{- 
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-''450iv) 


exp{~ 




exp{ 


-T450w) 


exp( 


-TS50w) 
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-T450iv) 
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SBssoN ea;p(— rssoiv) 

Given that we have shown we are coupling to sources larger in scale 
than 60 arcseconds, we could further correct these values by the 
ratio 

V450(cloud) %50(60") 



VSbOicloud) ^7450(60") 



(B13) 



where i] is the coupling. However, given uncertainties in the cloud's 
angular size (principally due to uncertainties in its height), and the 
fact that the ratio of the coupling values is likely to vary less signif- 
icantly than the individual values, we choose not to. 



APPENDIX C: REFRACTION 
CI Refraction Code at JCMT 

The primary purpose of the current JCMT refraction model is to 
provide coarse corrections to pointing based on local atmospheric 
conditions. It assumes a well-behaved, non-turbulent atmosphere. 
To first order, refraction displaces images towards the zenith by an 
amount 



R — A* tan(z) 



(CI) 



where z is the zenith distance, and yl is a function of the local 
atmospheric parameters : 

(i) T - Temperature (above some mean, chosen as 4 K for Mauna 
Kea) 



(ii) p - % Pressure change from the MK standard of 624 mb 

(iii) h - % Humidity 

Grids of atmospheric models were generated and refraction 
calculated at several zenith distances by integrating through the at- 
mosphere. Simple functional forms of the results were sought that 
would enable calculation of refraction from local atmospheric con- 
ditions. The optical- and millimetre- A-terms in the formula above 
take principally these forms (units are arcseconds) : 

A = 35.893 - 0.00067(;i - 20) - 0.135(r - 4) + 0.371p (C2) 

at 0.55 ii.m and: 

A = 36.800 + 0.0768(/i - 20) - 0.0294(r - 4) + 0.371p (C3) 
at 1 mm. 

- i.e. optical refraction is dependent mostly upon temperature, 
while submillimetre refraction is dependent mostly upon humidity. 
The 1 mm formulation used at JCMT contains additional, less sig- 
nificant terms. Its values were compared to the integration results 
throughout the grid and were found to be accurate to better than 
1 arcsec for zenith distances less than 80 degrees under all but the 
most extreme conditions. 



C2 Expectation Values of Signal Levels in the Presence of 
Refraction 

This appendix describes the use of the pointing rms values deter- 
mined from the SAO phase monitor to give an indication of the 
expectation values of the signal levels, relative to the values for no 
refraction 'noise'. The assumptions are 

(i) the refraction noise can be described by a Gaussian process, 
which in two dimensions gives the Rayleigh distribution: 



■ exp 



2cr2 



(C4) 



where Vg is the probability of the refraction having a value of 
arcsecs, and a is the rms refraction 'noise'. 

(ii) the telescope beam can also be defined as a Gaussian: 



S{e) = So exp 



-4;n(2)- 



|2 

jwhm 



(C5) 



where S is the signal level which would be measured for a move- 
ment of the telescope of 6 arcsecs, and Of^hm is the FWHM size 
of the beam in arcsecs. 

The expectation value of the signal level, relative to the no 
refraction noise case, is given by 



< 5" > = 



j^" S Ps dS 



where Pg is the probability distribution for the signal. 
Probabilities transform directly so that 

PsdS = Pe{-de) 



(C6) 



(C7) 



where use is made of the fact that the probability of S decreasing is 
a function of the probabilit y of 9 increasing. 
For simplicity, we rewrite (CSl) as 



S{e) = Soexp 



"2/32 



(C8) 
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where 13' = 4^ . 

^ Bln2 I — -I 

Now, from equation (|Cq), 

ds_ _ _e_ 

de ^ /32 

using this and equations (|c4|),(^8|), 

Ps 



(C9) 



P' 1 



S 



5^ 

So 



(CIO) 



It should be noted that, if 0' = a^, t he probability is uniform. 
If we now evaluate ( C6 ) using ( K ) we find that 

< S > 1 1 



So 



1 + 



J- + 32 



(CU) 



In the case of a perfect 15-metre telescope, one can write for Oj\„hn 
in arcsecs : 



^ fwhm 



[14.025 + 0.1856Tb]A 



(CI 2) 



where Te is the ill umination edge ta per in decibels and A is the 
wavelength in mm (Goldsmith 1987). The edge taper determines 



the level of ground radiation accepted by the feed, and at JCMT 
Te ~ 7.5 dB . In the case of a telescope with surface imperfections 
one must use the measured value of the full width half-maximum 



beam size in ( CI 1 ) 



If we assum e a 14.2" beam at 850 ^m and a 7" beam at 
450 [im, then ( CI 1 ) becomes 

1 



1 + 0.0275a2 
and 

1 



for the 14.2" beam 



for the 7" beam 



1 + 0.1132a2 
where a is the rms refraction noise in arcseconds. 



(C13) 



(C14) 



